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Abstract 



We have studied spin excitations in nanosized magnetic films in the Heisenberg model 
with magnetic dipole and exchange interactions by the spin operator diagram technique. 
Dispersion relations of spin waves in thin magnetic films (in two-dimensional magnetic 
monolayers and in two-layer magnetic films) and the spin-wave resonance spectrum in Al- 
layer structures are found. For thick magnetic films generalized Landau-Lifshitz equations 
are derived from first principles. Landau-Lifshitz equations have the integral (pseudodif- 
ferential) form, but not differential one. Spin excitations are determined by simultaneous 
solution of the Landau-Lifshitz equations and the equation for the magnetostatic potential. 
For normal magnetized ferromagnetic films the spin wave damping has been calculated 
in the one-loop approximation for a diagram expansion of the Green functions at low 
temperature. In thick magnetic films the magnetic dipole interaction makes a major con- 
tribution to the relaxation of long- wavelength spin waves. Thin films have a region of low 
relaxation of long-wavelength spin waves. In thin magnetic films four-spin-wave processes 
take place and the exchange interaction makes a major contribution to the damping. It is 
found that the damping of spin waves propagating in magnetic monolayer is proportional 
to the quadratic dependence on the temperature and is very low for spin waves with small 
wavevectors. Spin-wave devices on the base of nanosized magnetic films are proposed - 
tunable narrow-band spin-wave filters with high quality at the microwave frequency range 
and field-effect transistor (FET) structures contained nanosized magnetic films under the 
gate electrode. Spin-wave resonances in nanosized magnetic films can be used to construct 
FET structures operating in Gigahertz and Terahertz frequency bands. 

75.10.Jm; 75.30.Ds 

Heisenberg model, diagram technique, spin waves, nanosized magnetic films, relaxation, 
spin waves devices in Gigahertz and Terahertz frequency bands 

1 Introduction 



Nanosized magnetic films are of great interest due to their perspective applications in spin- wave 
devices. At present, the most important spin waves - microwave filters, delay lines, signal-to- 
noise enhancers, and optical signal processors have been realized on the base of magnetic films 
of microwave thickness [H [21 13]. Nanosized films give us opportunity to construct spin- wave 
devices of small sizes and to design devices with new functional properties. Recently new 
applications of spin waves have been proposed - spin- wave computing [HE], spin- wave filtering 
using width-modulated nanostrip waveguides [6], and transmission of electrical signals by spin- 
wave interconversion in an insulator garnet Y3Fe50i2 (YIG) film based on the spin-Hall effect 
[?]. Spin-wave logic elements have been done on the base of a Mach-Zehnder-type interferometer 
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[21 El [H] and can be realized on magnonic crystals [S] . Using nanosized magnetic films, we have 
probability to construct array of logic elements of small sizes. 

In order to design new spin-wave devices based on nanosized magnetic films, it is necessary 
to determine the dispersion relations and damping of spin excitations in nanosized films. In the 
phenomenological model with the magnetic dipole interaction (MDI) and the exchange inter- 
action fLU\ mi [121 113] the magnetization dynamics in thick magnetic films is described by the 
Landau-Lifshitz equations, which are differential with respect to spatial variables. The differ- 
ential form of equations is postulated. In this connection, the following question arises: is this 
form of Landau-Lifshitz equations correct for nanosized films? Determination of the dispersion 
relations depends on the answer of this question. In phenomenological models the spin-wave 
damping is described by relaxation terms in Gilbert, Landau-Lifshitz, or Bloch forms [I3]. Prop- 
erties of intrinsic relaxation processes are not taken into account in these terms and, therefore, 
the calculated spin- wave damping may be incorrect. The above-mentioned leads us to the main 
question of the paper: what are the dispersion relations and damping of spin waves in nanosized 
films and can they be derived from first principles? In order to answer this question, we develop 
the Heisenberg model with the MDI and the exchange interaction. In the framework of this 
model we consider spin excitations in nanosized films, relaxation of spin waves, and generalize 
Landau-Lifshitz equations. 

The above-mentioned problems have not yet been investigated comprehensively. One of the 
cause of these problems is the long-range action of the MDI. The spin-wave relaxation and the 
spin-wave dynamics become dependent on the dimensions and shapes of ferromagnetic samples. 
In order to analyze the Heisenberg model with the MDI and the exchange interaction we use the 
spin operator diagram technique [111 [151 IISI [13 |T8] . Advantages of the spin operator diagram 
technique are: the opportunity to calculate the spin wave damping at high temperatures and 
more exact relationships describing spin-wave scattering and excitations in comparison with 
methods based on diagram techniques for creation and annihilation magnon Bose operators 
|2ni[nil22l[2Sll21l2a[2ni[2Z]. In[ISll2B]the spin operator diagram technique is generalized for 
models with arbitrary internal Lie-group dynamics. 

In section 2 we consider spin operator diagram technique for the Heisenberg model with 
the MDI and the exchange interaction. Spin-wave excitations are determined by poles of the 
"P-matrix - the matrix of the effective Green functions and interaction lines. On the base of 
this diagram technique dispersion relations of spin waves in a normal magnetized monolayer 
and in a magnetized structure consisted of two monolayers and the spectrum of spin-wave 
resonances in a A^-layer structure are found (section 3). For thick magnetic films it is more 
convenient to present the P-matrix-pole equation describing spin-wave excitations in the form 
of the Landau-Lifshitz equations and the equation for the magnetostatic potential (section 4). 
Spin excitations are determined by simultaneous solution of these equations. Landau-Lifshitz 
equations are integral (pseudodifferential) equations, but not differential ones with respect to 
spatial variables. The reduction of Landau-Lifshitz equations to differential equations with 
exchange boundary conditions is incorrect and their solutions give dispersion relations differed 
from dispersion relations calculated on the base of integral (pseudodifferential) Landau-Lifshitz 
equations. In section 5 we consider spin-wave relaxation in thick and thin magnetic films. 
In thick films three-spin-wave processes take place and the MDI makes a major contribution 
to the relaxation of long-wavelength spin waves. Thin films have a region of low relaxation 
of long- wavelength spin waves. In this case, three-spin-wave processes are forbidden and the 
exchange interaction makes a major contribution to the relaxation process. Nanosized magnetic 
films with low relaxation spin waves are applicable to microwave spin wave devices. Tunable 
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narrow-band spin-wave filters with high quahty at the microwave frequency range and field- 
effect transistor (FET) structures contained nanosized magnetic films under the gate electrode 
are proposed in section 6. Spin- wave resonances of nanosized magnetic films can be used to 
construct FET structures operating in Gigahertz and Terahertz frequency bands. 

2 Heisenberg model with magnetic dipole and exchange 
interactions 

2.1 Spin operator diagram technique 

Let us consider the Heisenberg model with the exchange interaction and the MDI on a crystal 
lattice [TTIIIH]. The exchange interaction is short-ranged and the MDI is long-ranged. Operators 
= ± iS^, satisfy the commutation relation 

[S^{l),S^{n = S^l)6rr 
[S%l),S-{l')] = -S~il)6^,^, 

[S^l),S-il')] = 2S^il)5^^„ 

where 1 = ri, 1' = rl ' is the abridged notation of crystal lattice sites. 
The Hamiltonian of the Heisenberg model is 

1 1 1,1' 

where H {H || Oz) is the external magnetic field, is the auxiliary infinitesimal magnetic 
field, /i = — , +, z. It is supposed that the summation in ([T]) and in the all following relations is 
performed over all repeating indices n, v. The summation is carried out over the crystal lattice 
sites 1, 1' in the volume V of the ferromagnetic sample, g and ixb are the Lande factor and the 
Bohr magneton, respectively. J^uiJ — 1') = JufiiJ' — 1) is the interaction between spins, which 
is the sum of the exchange interaction 7^,^ and the MDI 

7^,(1-1') = ^(1- 1') - An{gf^ByV,^f- ^^1^^^,,^^, , (2) 
where $(r — r ') is determined by the equation 

A$(f-f') = 5(r -f'), (3) 

In the 3-dimensional space $(r — r ') = —1/An\f — f'\ and the MDI term in the Hamiltonian 
([1]) can be written as 
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For the following calculations of spin-wave dispersion relations in magnetic films we use more 
convenient form of the MDI determined by relations ([2]), ([3]). 

Spin excitations, interaction of spin waves, spin wave relaxation and other parameters of 
excitations in the canonical spin ensemble are determined by the generating functional [T4t [29| 

Z[h] = Spexpl- I3n{h)] 
E E /■ ■ ■/ . . . , n, Ti, . . . , r„) V(r, Ti) . . . r„) dr^... dr^, (4) 



n=0 1 







where P = 1/kT, k is the Boltzmann constant, T is the temperature, h = {hf^.}. Coefficients 
Qiii,...,iin gj-Q proportional to the temperature Green function without vacuum loops 

G^^-^"il . . . , n, n, . . . , = {{TS^^ (1, n) . . . S^" (n, r„))) 



(5) 



where S"'{n, r) = exp{T'H)S°'{n) exp(— rH) are the spin operators in the Euclidean Heisenberg 
representation, r G [0, /?]. T is the r-time ordering operator. Variable r is added in the auxiliary 
field in order to take into account T-ordering. ((...)) denotes averaging of spin operators 
calculated with exp(—/37/)/Spexp(— The symbol Sp denotes the trace. 

The frequency representation of the expansion (jl]) is more convenient for calculations. 
The Fourier transforms of Q/^i'---'/^" are defined in terms of the Matsubara frequencies co^^^ 
2Trmi/hf3, . . ., = 27nnn/hf3 [30] (mi, . . ., m„ are integers) 



mi 
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= J...J Q^^'-'^"(r, . . . , n, n, . . . , r„) exp[-ih{ul^,lT, + ...+ a;gr„)] dn . . . dr^,. (6) 



The coefficients Q^^^^■■■'^^" can be expanded with respect to the interaction J^j,(l — 1') ([2]) [T71 
[TS| [m [151 IISI EH]- Each term of this expansion is represented by a diagram constructed of 
propagators, vertices, blocks and interaction lines. 

1. Propagators. Spin propagators 

D±(r,l',u;^) = ^%— , (7) 

where po = Pqij^bH, are determined for the spin ensemble without any interaction between 
spins. The propagators D±{1, 1', Um) are represented by directed lines in diagrams (figure [T](a)). 
The directions of arrows show the direction of growth of the frequency variable Um- 

2. Vertices. There are five types of vertices (figure [T](b)). Vertices a, b are the start and 
end points of propagators, respectively. In analytical expressions of diagrams the vertex a 
corresponds with the factor 2 and the vertex b with the factor 1. The vertex c ties three 
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a. D,(1,1',(dJ 




Figure 1: (a) Propagators D±, (b) vertices, (c) block with isolated parts and (d) interaction 
lines V;(°). 



1' 
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propagators and corresponds with the factor (-1) in analytical expressions. The vertex d with 
the factor 1 is defined as a single vertex. The vertex e ties two propagators. The factor of the 
e- vertex is equal to (-1). 

3. Blocks. Blocks contain propagators and isolated vertices d (figure [T](c)). Propagators can 
be connected through vertices c, e. In analytical expressions of the diagram expansion each 
block corresponds with the block factor B^'^~^\pq)^ where k, is the number of isolated parts in 
the block. The factor B^'^~^\pq) is expressed by partial derivatives of the Brillouin function Bg 
for the spin S with respect to po 

B{p,) = {{S^))o = SBsiSpo) 

where ((. . .))o denotes the statistical averaging performed over the states described by the 
Hamiltonian H ([1]) without the interaction J^^^, between spins. Bs{x) = (1 + 1/2S') coth[(l + 
l/2S)x]- {l/2S)cot\i{x/2S). 

4. Interaction lines. The interaction line V^^°^(l — 1', cjm) = PJiiui^—^') connects two vertices 
in a diagram (figure [T](d)). The correspondence between the first index fi of the interaction line 
V^^^^ and the vertex type is the following. (1) If /i = — , then the left end point of vl^J is bound 
to the vertex a; (2) if /i = +, then this end point is bound to the vertices b or c; (3) if fi = z, 
then the end is bound to the vertices d or e. The analogous correspondence is satisfied for the 
rig ht end z/of 

Coefficients Q'^i'---''^" in the expansion (j3]) in the frequency representation ([6]) are the sum of 
N topologically nontrivial diagrams J2n Qtv ' The general form of the analytical expression 
of the diagram in the frequency representation is written as [T71 [TSl HU [T51 [IS] 

Q%^'-'''"{1, . . . ,n,a;S, • • • = (-l)"2-'"^n5"''"''(Po) ft % 



X E EK^;^^-l^c^n^Jx...xv;(°)(^-^^c^^J^D4s,.-^c^^J^4^^^^^^ 

r....fe mi g^g' V \r&v J 

i'...k' 

where 1, . . . , n, ojI^^ , • • • , ojml the external lattice and frequency variables corresponded to 
the auxiliary fields /i^. in the expansion (jl]). is the number of a- vertices in a diagram. L is 
the number of c and e-vertices. Pk is the number of topological equivalent diagrams. 2k is the 
number of vertices connected with k interaction lines V^^°^ . . . V^^^ . The product Hi is performed 
over all blocks of a diagram, kj is the number of isolated parts in block /. The term Yl'^'^ , 
denotes that all isolated parts in block / are determined on a single crystal lattice site. Id is the 
number of propagators in a diagram. J„ is the number of vertices in a diagram. J2m, denotes the 
summation performed over all inner frequency variables. The term Yii'' S {J2rev (^^^mr) gives 
the frequency conservation in each vertex v, i.e. the sum of frequencies of propagators and 
interaction lines, which come in and go out from the vertex v, is equal to 0. The vertex d can 
be connected with the single interaction line. In the analytical expression this corresponds to 
the factor 6{PhuJm)- The lattice variables s, s' of propagators D_ can be inner or external. In 
the first case, end points of propagators are connected with the end points {1, 1', . . . ,k, k'} of 
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interaction lines . . . V^;^^ and the summation Z^m, is performed. In the second case, 

i'...k' 

end points of propagators are not connected with interaction hnes. 

The first approximation of the diagram expansion (jl]) is the self-consistent field approxi- 
mation, in which the effective field acting on spins is derived and the self-consistent field H^^^ 
induced by the neighboring spins is taken into account [m [171 [IE]. This leads to the substi- 
tution Pq ^ p = l3giJiBHf' in the propagator in relation ([7]). The self-consistent field is the 
sum of exchange and magnetic dipole self-consistent fields, H^f> = Hjf^^^^ + Hj^^^ , where 



^f'^'Hi) = (^7/ii.)-'EV(i-i')((^'^(iO)) 

1' 

(10) 

r' = l' 

The second approximation of the expansion (jl]) is the approximation of the effective Green 
functions and interactions. In this approximation, the poles of the matrix of the effective Green 
functions and interactions are determined and the dispersion curves are obtained. The next 
terms in the diagram expansion determine the imaginary and real corrections to the poles of 
the matrix of the effective Green functions and interactions. The imaginary parts of the poles 
give the relaxation parameters of spin excitations and the real parts determine the corrections 
to the dispersion curves. In the next section we consider the approximation of the effective 
Green functions and interactions. 



1' 



2.2 Effective Green functions and interaction lines 



In the framework of this approximation the matrix of the effective Green functions and effective 
interactions V = || Pas (1, 1', i^m)|| is introduced [171 IlH]- We compose the P-matrix from 
analytical expressions of connected diagrams with two external sites. These sites are end points 
of propagators, single vertices d, or end points of interaction lines. Accordingly, multiindices 
A = (afi), B = (6z/) are the double indices, where = { — , +, z} and indices a, b point out 
that A, B belong to a propagator or to a c?- vertex (a, 6 = 1), or belong to an interaction line 
(a, 6 = 2). The zero-order approximation V^^^ of the P-matrix is determined by the matrix of 
the bare interaction V*^°^ = ||V^^°Hl — l',Wm)|| and by the two-site Green functions ^ in the 
self-consistent-field approximation Q'^'^^ = given on a crystal lattice site 



V 



(0) 



(0) 



\P, 



(0) 



(1m)(2i/) I 



|p(0) II 



V ; wv^Jw) 



where 



i 















2B(p)D+(r,r',a;, 




25(p)D_(r,r',a;„ 











(11) 
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with the propagator ([7]), in which the substitution Po ^ p = PqubH^^^ is performed. 

The P-matrix is obtained by means of the summation of the 7^*^°^ -matrix - the summation of 
all diagram chains consisted of the bare Green functions G^^} and the bare interaction lines V^^°^ 
(figure E]). These chains of propagators and interaction lines do not have any loop insertion. 
Analytical expressions of the considered diagrams can be written in accordance with relation ([9]). 
The summation gives equation of the Dyson type, which forms the relationship between pC^)- 
and P-matrices 

p = p(o) +p(^p(o)^ (12) 

where 

\ £ J 

£ = W^fiuW is the diagonal matrix. 

The "P-matrix consists of the two-site effective Green functions Q = \\G^^\\ = Q^^\£ — 
V(o)^(o))-i, where = effective interactions V = ||V;,^|| = V^^^^ - gio)y{o))-\ 

where V^,y = P[2^i){2u), and intersecting terms P(i^)(2i.), -P(2M)(ii/) (figure E]). The effective Green 
functions, effective interactions and intersecting terms are denoted in diagrams by directed thick 
lines, empty lines and compositions of the thick line - empty line, respectively. The P-matrix 
determines the spectrum of quasi-particle excitations in the spin ensemble. Spectrum relations 
for spin excitations are given by the "P-matrix poles - by zero eigenvalues of the operator 
1 — crP^°) or, equivalently, by £^ — V'-^-'^*-'^'' under the analytical continuation 

iujm OJ + ie sign u 

8{l3huJm) = ^ ^7——^. T (^^+0). (13) 

pri[uj + le signcjj 

Since, zero eigenvalues of the operator 8 — V*-"-*^ can corresponds to different eigenfunc- 
tions and can determine different excitation modes, we introduce the spectral parameter A for 
eigenfunctions h'^^\l,0Jm) of the operator 8 — V*^°-'^*^°-'. The spectral parameter A can be dis- 
crete or continuous. Taking into account the above-mentioned, we get the equation describing 
spin-wave excitations 



IM'' 



iijJrn — >w+ iesigntj 



0. 



(14) 



3 Spin waves in nanosized magnetic films 
3.1 Spin- wave equations for magnetic films 

Let us consider spin waves with the wavevector g* in a normal magnetized film consisted of N 
monolayers at low temperature, x-, y-ax.es are in the monolayer plane and the 2;- axis is normal 
to monolayers. The external magnetic field H is normal to monolayers and is parallel to the 
^-axis. At low temperature derivatives of the Brillouin function in B^'^\p) in relation ([HD tend to 
exponentially with decreasing temperature. Thus, it follows that diagrams containing blocks 
with isolated parts can be dropped, the Green function G^^) in relation (fTTj) is negligible and 
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Figure 2: (a) Definition of the effective Green functions P(i^){iu) = G^u via the bare two-site 
Green functions G^^}. (b) Definition of effective interaction fines P(2^){2p) = V^u- (c) Definition 
of intersecting terms P(i^)(2i/), P{2^l)(lu)■ 
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only the Green functions G_\_, G+_ are taken into account in equation f|T4|) . Indices /i, u of 
interactions V^^^^ in equation f lT^ are { — ,+}• We suppose that on monolayers spins are placed 
on quadratic crystal lattice sites with the lattice constant a and spin orientation is parallel to the 
z-axis. The exchange interaction acts between neighboring spins and is isotropic between spins 
in monolayers, 2/_+ = 2/+_ = I^z = Iq, and between neighboring layers, 2J_+ = 2/+_ = /^^ = 
Id- Then, the Fourier transform of the exchange interaction with respect to the longitudinal 
lattice variables Ixy is 

/(g,i,-i;)= ^ j(r^^_r;^,i^_i;)expHg(r,,,-rg] 



= /(0, 1^ - 1^) + 2Jo[cos(g^.a) + cos{qya)]5i^r^, 

where l^y, I'^y are crystal lattice sites in a monolayer, 1^, 1'^ are 2;-positions of layers, q= {q^, qy) 
is the longitudinal wavevector in monolayers, J(0, 1^ — 1^) is the exchange interaction at g'= 0, 
which is equal to Id between spins of neighboring layers. The corresponding exchange part of 
the interaction line V^^) = V^^^'^'^) + V^^fp'^ (figure [^d)) is 

v^r^\qA.-i':) = pmu-i'.). as) 

The MDI part V^^^*^^ is determined by the Fourier transform of equation ([3]) 

{ -<f + ^ ^i^^ z-z') = S-H{z - z') 



dz 



with the solution 



m u - I'j = m z - ^')L=i.,.'=i; = ^^^p(-^ii^ - I'-D' (16) 

where Sa = a? , q = \q\- According to the solution f|T6|) . the corresponding MDI part of the 
interaction line is 

V!^^^\q. U - 11) = -^"^(^g^)'^^^- exp(-g|l. - i;|), (;., . = {-, +}). (17) 

where 

Taking into account relations f|T5l) and f|T7|) . from equation f|T4l) we obtain equations for spin- 
wave modes with the wavevector q in A^-layer magnetic films 

h^^\q, lz,u^) - J2 [v!l'J{q, 1. - 1Lu;^)GL°1(1;, l',,oojh':^\q, 
+ V;(?(g,l.-l',,a;JGyi(i;,l',,u;J/iL')(g,i;,u;jl , _„,, = 0, (18) 

where 

^(0) _ 2E(p),5i,i. 

(+-) p ± ipnUrr, 
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Figure 3: Dispersion curve of spin waves propagating in the normal magnetized monolayer film 
with cubic lattice (a = 0.4 nm) at the sum of magnetic fields H + if (™) = 3 kOe. Exchange 
interaction Iq is 0.085 eV. 

A = 1, . . . , iV is the mode number, V^^iq, U - 1'., i^m) = VjiT'^Xq, U - I'J + V;(f^')(g', 1, - 1^), 
IX, u = { — ,+}. Eigenvalues of equations (ITS!) give dispersion relations of spin waves. In next 
sections we find spin-wave dispersion relations for the cases of monolayer and two-layer films 
and spin-wave resonance relations for the case of A^-layer structures. 

3.2 Spin waves in magnetic monolayer 

Dispersion relations of spin waves in normal magnetized monolayer are determined by the 
determinant of equations ( |T8l) for variables h^l\ Iv^K Taking into account relations ( |T5l) and 
(Tfl). we find 

uj\q) = n{q}[n{q} + 27r7a^g], (19) 

where 

n{q) = 7(iJ + iJ^'")) + ^^M^[2 - cos(g,a) - cos(g,a)], 

7 = qiib/T^ is the gyromagnetic ratio, if^™) is the depolarizing magnetic field ( TTOl) . cxm = 
qUbB^p) / Sa is the surface magnetic moment density, q = (ql + qyY^"^- As one can see from 
relation ( !T9|) . in monolayer films spin waves have the one-mode character. Figure |3] presents the 
dispersion curve (fT9|) of spin waves propagating in the monolayer film with the lattice constant 
a = 0.4 nm. The spin-wave wavevector q is parallel to the x-axis {qx = q, qy = 0) and is in the 
range [0,7r/a]. Calculations have been done for the exchange interaction between neighboring 
spins Iq = 0.085 eV, B{p) = 1/2 at the sum of magnetic fields H + H^"^^ = 3 kOe. The 
exchange interaction makes a major contribution to the dispersion. The relatively weak MDI 
is significant for the dispersion at small values of the wavevector q. 
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Figure 4: Dispersion curve of spin waves propagating in the normal magnetized two-layer 
magnetic film with quadratic lattice (a = 0.4 nm) at the sum of magnetic fields H + H^"^^ = 
3 kOe. Exchange interactions Iq = Id = 0.085 eV. Distance between monolayers d is equal to 
the lattice constant a. 1, 2 are the first and the second modes of spin waves, respectively. 

3.3 Spin waves in two-layer magnetic film 

Let us consider spin waves in a normal magnetized structure consisted of two monolayers of 
the quadratic lattice with the lattice constant a. The distance between layers is equal to d and 
the exchange interaction between spins of layers is Id- Dispersion relations are determined by 
eigenvalues of equations ( !T8l) for variables , h^^\ h^^\ fv^ and can be written as 



±2 

where 



h 

2B{p)Id 



27r7(Tmgexp(-grf) 



(20) 



m = ^{H + + ^^M^[2 - cos(g.a) - cos(g,a)] + 

n = 1,2 is the mode number, q = (q^ + qyY^'^. For the first mode spins in different layers 
change their orientations in-phase. In this case, spin waves of the first mode correspond to 
spin waves in monolayer f[TI?|) . For the second mode spins in different layers change orientations 
in-anti-phase and the energy of the spin wave with the given longitudinal wavevector q is higher 
than the energy of the spin wave of the first mode. Dispersion curves of spin waves determined 
by relations (l20l) are shown in figure HI Spin waves propagate along the x-axis. Calculations 
have been done for the exchange interactions Iq = Id = 0.085 eV and for the distance between 
layers d = a = 0.4 nm at the sum of magnetic fields H + = 3 kOe. 
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3.4 Spin- wave resonance in A^-layer structure 



In this section we consider spin- wave resonance in a structure consisted of uniform monolayers 
with the exchange interaction Id between spins of layers and with the distance d between layers. 
Spin-wave resonance is the limit case of a spin wave, when the longitudinal wavevector g — )■ 0. 
Therefore, the MDI terms V^^'^\q, 1^ — 1^) in equations ( |T8i) can be dropped, equations with 

variables h^^^ and /i^^ are separated and eigenvalues are determined by zero values of the 
determinant (we write the determinant for equations with the 



X det 



-y(o)(21) 




-yW(12) 





(G(o)-i(2) -r(°)(22)) -y(o)(23) 

-\/W(32) (G'W-H3)- V(°)(33)) 



v 



where V'^^\ij), G^^\i) are the abridged notation of v}!^^^^\q,iz — jz,^m)\q=o and G^^\_ 
at ium uj + iesignu, respectively. are indices of layers. Taking into account that spins 

of outer layers {i = 1, N) interact with spins of one inner layer and spins of inner layers interact 
with spins of two layers and introducing the variable for inner layers in the determinant D*^"*") 



G(o)- 



X 



V(0)( 



It 



h 



.[CO - ^{H + H^^^)] - 2 it^l,N, j=2±l). 



we obtain that the spin-wave resonance spectrum is determined by roots of the polynomial 



Rn{x) 



(x + 1) 
1 
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X 1 
1 (x + 1) 



= {x + lfPN-2{x) - 2{x + l)Piv-3(a;) + Pn-a{x) = {N > 2), 

where P-2{x) = — 1, P-i(x) = 0, Pq{x) = 1, Pn{x) = xPn-i{x) — Pn-2{x). Polynomial Rn{x) 
has roots 

a:(") = -2cos(^), 

where ?t, = 0, 1,...,A^ — 1. Taking into account the form of the roots x^'^\ we can introduce the 
transverse wavevector g^*^^ = nn/Nd. Then, the spin- wave resonance spectrum can be written 
as 



(if + ijM) + ^^|H:^[l_cos(gi")rf)]. 



(21) 
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Figure 5: Spin-wave resonance spectrum w^"-' (n = 0, 1, . . . , A^ — 1) for the structure with = 40 
layers, g^"^ is the transverse wavevector, d is the distance between layers, Id is the exchange 
interaction between spins of layers. 

For the first mode {n = 0) spins in different layers change their orientations in-phase. For the 
highest mode {n = N — 1) spins in different layers change orientations in-anti-phase and the 
energy of spin-wave resonance is highest. Figure [5] presents the spin-wave resonance spectrum 
for the structure with = 40 layers. One can see that at low values of the transverse wavevector 
the resonance spectrum is proportional to the quadratic dependence on q[^\ 

4 Landau-Lifshitz equations and spin-wave excitations 
in thick magnetic films 

4.1 Linearized Landau-Lifshitz equations 

Equations ([HI), (ITSl) describe spin- wave excitations. Solutions of these equations for magnetic 
samples of great volumes and for thick A^-layer magnetic films with ^ 1 become difficult, be- 
cause determinants of equations ( |T4l) . (JTSll have high orders. In order to overcome the difficulty 
and to find spin- wave spectrum for these samples, we derive Landau-Lifshitz equations p/Tj il8j . 
Dispersion relations for spin excitations are determined by the P-matrix poles f|T2|) which coin- 
cide with poles of the matrix Q of effective propagators. Accordingly, the dispersion relations 
can be derived from the eigenvalues of equation 

g = ^(0) + ^(V^''^'^) + v('^'P^)g(^\ (22) 

where Q^^^ = \\G^^-}\\ is the matrix of bare propagators (ITTl) . Since the considered interaction 
is the sum of exchange and magnetic dipole interactions, we can obtain the eigenvalues and 
eigenf unctions of equation fl22l) by a two-step procedure. In the first stage, we perform the 
summation of diagrams, taking into account the exchange interaction, and find the propagator 
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matrix ^(i) = ||GW|| 

^(1) = ^(0) ^^(0)y(exch)^(l)_ ^23) 

In the second stage, the summation of diagrams with dipole interaction hnes is performed. 
This gives the equation for the matrix Q of effective propagators expressed in terms of the 
matrix Q^^^ 

g = + gv^'i^p)gW. (24) 

Thus, the solution of equation ( l22l) . which determines the matrix g, is equivalent to the 
solution of equations (!23|) . (^^. After the performed two-step summation, equation (fT4|) for 
eigenfunctions h^^^ is written in the more convenient form 



UJr^ 



= 0. (25) 



The solution of simultaneous equations f l23|) . f l25|) gives the dispersion relations for spin 
excitations. These equations can be reduced to linearized Landau-Lifshitz equations in the 
generalized form and the equation for the magnetostatic potential. In order to perform this 
transformation one needs to make a transition to the retarded Green functions. We transform 
matrix equation (123|) to equations describing small variations of the magnetic moment density 
(or the variable magnetization), m^. The variable magnetization rrii, under the action of the 
magnetic field h^, which is generated by the MDI V^^^^\ is given by the retarded Green functions, 
which are determined by the analytical continued values of the propagator matrix g^^^ 

m.(r,c.) = ^^EG2(r,f',c.„ 



a 



h,{l',u), (26) 



where Va is the atomic volume. The analytical continuation ium OJ — is defines the retarded 
Green functions. /ip(l,a;) is the field of the magnetic dipole-dipole interaction acting on spins. 
By multiplying matrix equation fl23|) by from the left and by hp from the right, performing 

the analytical continuation iujm u — ie, 6{f3hujm) [(3h{uj — ie)]^^ and taking into account 
relation f l26|) . we get matrix equation f l23|) in the form of simultaneous equations 

E[G(°)-^(r, l',u;) - /3V(r- f')K(f',a;) = ^^^^hpHu). (27) 

For isotropic exchange interaction, 2/_+ = 2/_|__ = I^z = I, equations (1271) have the form 

E±m±{l, to) = 2-fM{l)h^{l, to) (28) 



E.m^itu) = ?^^M{T)K{l,u), (29) 
B{p) 

where M(l) = g^BB{p) /Va is the magnetic moment density at the low-temperature approxi- 
mation. We say that the operators E±, E^'- 
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+ ^ E / [^0) - m] exp[^g(l - V)]m^{V, oj) d'q 

^ 1*' Vt, 

are Landau-Lifshitz operators. I{q) = X]r-^(1) exp(— igl) is the Fourier transform of the ex- 
change interaction with respect to the lattice variables. The field H^"^\l) is defined by relation 
( 1T0|) and depends on the magnetic moment density M(l); Vb = (27r)'^/V^ is the volume of the 
first Brillouin zone. Equations (1281) . fl29l) have the generalized form of the Landau-Lifshitz 
equations [13]. Solutions m± of equations fl28|) depend on temperature, because /3 = 1/kT 
is contained in the variable p of the function B{p) ([8]), through which the magnetic moment 
density M(l) is expressed. Equation ( l29l) describes longitudinal variations of the variable mag- 
netization under the influence of the field hz- At low temperature the derivative of the function 
B^^^(p) tends to and the longitudinal variable magnetization is negligible. 

From the form of the magnetic dipole interaction in relations ([2]), ([3]) it follows that the 
field hi, in relation (126|) is magnetostatic, i.e. it is expressed in terms of the magnetostatic 
potential ip: = — V^y?. We transform equation (125]) to the equation for the magnetostatic 
potential y)(r, cj). Taking into account relation (1261) and the explicit form of the magnetic 
dipole interaction in relations ([2]), ([3]), performing the derivation V^, the analytical continuation 
iiOm u: — is and the summation of equation fl25|l over the index /i, we obtain the equation 
expressed in terms of rriy 

- A(/?(f,u;) +47rV^m^(l,a;)|j'^,.= 0. (30) 



Thus, in consideration of the Landau-Lifshitz equations (1281) . (1291) . the dispersion relations of 
spin excitations are given by eigenvalues of equation fl30l) . 

4.2 Exchange boundary conditions 

If the scale of the spatial distribution of the variable magnetization m^(l,u;) and the sample 
size are much greater than the lattice constant a, then the sum over the lattice variables X^f 
in Landau-Lifshitz operators E±^ Ez can be converted into an integral over the sample volume 
J d^r. Let us consider the case when the temperature is low and the Fourier transform of 
the exchange interaction is I{q) = 7(0) —wq^. Then, we obtain that — )■ and equation fl29l) 
is dropped. The operators E± are pseudodifferential operators of order 2 

E±m±(f,u) = [-f{H{f) + H^"'\r)) ± uj]m±{f,uj) 
47T'~faM{f) 



(2n 



l3 



q exp[iq{f—r')]m±{r',uj)d qd r, (31) 

V Vi 

where a = wVa/4:7T{gnBy is the exchange interaction constant, V is the volume of the ferro- 
magnetic sample. In [TOl IHl |13l |33l |3ll |35l |36l EI] the pseudodifferential Landau-Lifshitz 
operators are reduced to the differential operators with respect to spatial variables 

E±{f,u) = -i[H{r) + H^'^\r) - 47raM(f)A] ± uo. (32) 
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For solvability of equations f l28p with differential Landau-Lifsliitz operators (152]) the ex- 
change boundary conditions are imposed 

on 

where ft is the inward normal to the boundary dV, and ^ is the pinning parameter. This reduc- 
tion to differential Landau-Lifshitz operators is not correct. Figure [6] presents exact dispersion 
relations of spin excitations given by eigenvalues of equations fl28|) . fl30|) with pseudodifferential 
and differential Landau-Lifshitz operators for the case of a normal magnetized homogeneous 
film with the thickness D. The dispersion relations of spin waves have the form 

Jn)2^^ = + n^.gV^o"^'), (33) 

where n = 1,2,3,... is the mode number, q = {qx,(ly) is the two-dimensional longitudinal 
wavevector of spin waves, ^=1^1, = 7(if — ArcM + ArcaMq^"^'^), ^Im = 47r7M, gg"-* = 
(g^ + qi^^'^Y^'^, qi""^ is the transverse vector. The magnetostatic potential over thickness z G 
[—D/2, D/2] of the magnetic film is 

= (27r)-V(")-^/2gxp(^g,.x + ^g,i/)cos[gl")^ + 7r(n- l)/2], (34) 

where = D/2 + g/gS"^^ 

For the case of pseudodifferential Landau-Lifshitz operators flHTl) . the transverse wavevector 
g|"^ is closely connected to the longitudinal wavevector q by the relation 

2cotg(")D = ^-^. (35) 
? qz 



For the case of differential Landau-Lifshitz operators (|32|) . the transverse wavevector is 
determined by the exchange boundary conditions and is given by the equation [TOj [13] 

2cotg(")D = % ^. (36) 

s qz 

Dispersion relations ( 133|) of the first spin-wave mode propagating in the YIG film of the 
thickness D = 0.5 /zm with AttM = 1750 Oe, a = 3.2-10"^^ cm^ at the applied magnetic 
field H = 3000 Oe are shown in figure [6] for the transverse wave vector q^^ (|35l) and for the 
transverse wave vector q[^^ ( l36l) with different pinning parameters ^. One can see that there 
does not exist any pinning parameter ^, at which the curve A calculated on the base of relation 
( 135|) coincides with the curves calculated on the base of the exchange boundary conditions. 
Thus, we conclude that the reduction to differential Landau-Lifshitz operators and the use of 
the exchange boundary conditions are incorrect. 



5 Spin- wave relaxation 

In this section we answer the question: what is the value of spin-wave relaxation in the model 
with magnetic dipole and exchange interactions derived from first principles? The answer 
depends on the ratio of the spin-wave energy to intervals between modes of the spin-wave 
spectrum and is different for thick and for thin magnetic films. In thick films the spin-wave 
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Figure 6: Dispersion curves of the first spin- wave mode propagating in the YIG film of the 
thickness D = 0.5 fim with AtcM = 1750 Oe, a = 3.2-10^^^ cm^ at the apphed magnetic 
field H = 3000 Oe. The curve A is calculated on the base of relation fl5^ for the case of 
pseudodifferential Landau-Lifshitz operators (!3T|) . Curves 1-4 are calculated for the case of 
differential Landau-Lifshitz operators ( 132|) on the base of relation (!36|) with different pinning 
parameters ^. (1) = 0.01, (2) 0.1, (3) 1, (4) 10. 



energy is greater than energy gaps between modes and a three-spin-wave process takes place. If 
the exchange interaction is isotropic, it cannot induce three-magnon processes and, therefore, 
the MDI makes a major contribution to the relaxation. We consider the spin-wave damping in 
thick films in the one-loop approximation. In thin magnetic films (for example, in nanosized 
films) the energy of long-wavelength spin waves is less than energy gaps between modes and 
three-spin-wave processes are forbidden. In this case, four-spin-wave processes take place, the 
exchange interaction makes a major contribution to the relaxation, and the spin-wave damping 
has lower values in comparison with the damping in thick films. We calculate the spin-wave 
relaxation for four-spin-wave processes in thin films for long-wavelength spin waves in the two- 
loop approximation. 



5.1 Spin- wave relaxation in thick films 

The spin-wave relaxation induced by a three-spin-wave process in normal magnetized homoge- 
neous ferromagnetic films is considered in [T71 [18] in the one-loop approximation for spin waves 
with small longitudinal wavevectors at low temperature. The relaxation is determined by self- 
energy diagram insertions to the P-matrix given by relation f ll2p (figure [7]). Damping 
of the j-mode excitation is defined by the imaginary part of the pole of the effective Green 
functions = P(i_)(i+) with insertions under the analytical continuation f|T3|) 
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Figure 7: Self-energy diagrams in the one-loop approximation at low temperature. B is deter- 
mined by relation (|8]). 



+ 



8B{p) 
where 



Pii-)(2z){i,qi,UJn)P{2z)ii+){k,q- qi,ujm - UJn)]N'^{j,q;i,qi; k,q- gl) (fqi 



iojm -^uj+ie sign ui 

(37) 



P(i_)(i+)(j, q, = 2plf (^^'^ + '^V^-l + ^^m) 



-P(l-)(22)(j!'f, l^n 



i2z){ 



pii) 



(, ,{j)2 I , ,2 \-l 



P 



'/Mi/ (j)2 



+ ,2) 



q± = ^{qxTiqy), 

is the Fourier transform of the exchange interaction, 



n 



E 

o"l,(T2,a-3 



sm 



Niji,qi;j2,q2; 33,(13) 



exp uXl^fc^Ofc ~ 



k=l 



is the block factor in the representation of the functions (134|) . /^■'•' = D/2 + q/q^ ^"^, (Tk = ±1; 
Ecri,<T2,o-3 denotes the summation over all sets {cti, (T2, cjs}. The spin-wave frequency u^^^ and 
the transverse wavevector q^/'' are determined by relations (1331) and (1351) . respectively. The 
damping A^-'-' increases directly proportionally to the temperature. 

Relation (1371) describes relaxation of the spin-wave j-mode caused by inelastic scattering on 
thermal excited spin wave modes. Relaxation occurs through the confluence of the j-mode with 
the fc-mode to form the i-mode. From the explicit form of the block factor in relation (1371) it 
follows that the confluence processes take place when the sum of mode numbers j + i + k is equal 
to an odd number. The confluence processes are induced by the MDI and are accompanied by 
transitions between thermal excited i- and fc-modes. Transitions take place when equation 
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Figure 8: Spin wave damping A^^-* = ^cj'^^Y'^^^'' of the first mode in normal magnetized YIG 
film with the magnetization 4nM = 1750 Oe and the exchange interaction constant a = 
3.2-10-12 cm^ at H = 3000 Oe, T = 300 K at different film thickness D. (1) D = 500 nm, (2) 
300 nm, (3) 200 nm, (4) 120 nm. 



= u;«(g W) - - g W) (38) 

has at least one solution g for the given g, i, j,k. Existence of solutions q^'^^ of equation f p8|) 
depends on the thickness of the magnetic film. With decreasing film thickness D, the density 
of dispersion curves of modes on the plane (cu, q) decreases and the frequency of the spacings 
between curves increase. The least frequency spacing occurs between the first [i = 1) and the 
third {k = 3) modes. Figure M shows the damping A^^^^ of the first spin wave mode versus 
the longitudinal wave vector q normalized by the film thickness D at different film thicknesses. 
Calculations have been done for a YIG film with the magnetization AttM = 1750 Oe and the 
exchange interaction constant a = 3.2- 10"^^ cm^ at if = 3000 Oe and T = 300 K. One can see 
that for the YIG film with the thickness D = 120 nm in the region qD < 0.14 the damping 
A^^^ is equal to due to the absence of transitions between modes. Thus, in thin magnetic 
films a low spin-wave relaxation region takes place. For the given j-mode this region appears, 
when the excitation frequency uj^^\q) is less than the difference uj^^\q^^^) — uj^^\q — g^*-*) at 
any values of the wavevector g**^**^. For the first mode co^^^ in the YIG film the low spin- wave 
relaxation region is shown in figure |9] at g — t- 0. If the film thickness D < D q{ijJ^^^ /2tx)^ then 
oj'^^\Qi) < uj^^\q'^^'>) — uj'^'^\q^^^) and the first mode has low values of the spin wave damping 
A«. 

5.2 Relaxation in thin magnetic films 

What is the value of spin wave damping in the low relaxation region in thin magnetic films? 
We consider four-spin-wave processes in the normal magnetized monolayer of the quadratic 
lattice with the lattice constant a at small longitudinal wavevector values q = {qx, qy) at low 
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Figure 9: Film thickness of YIG film versus the excitation frequency UJ^^\q)/2T[ of the first 
mode at the wavevector g* — t- 0. Low relaxation region of the first spin-wave mode exists for 
YIG films with the thickness D < Do{uj^^^ /27t). 

temperature. As isotropy of the exchange interaction can not forbid four-spin-wave processes 
and the value of the exchange interaction much greater than the MDI, only the exchange 
interaction will be taken into account in diagrams. We suppose that the exchange interaction 
acts between neighboring spins and is equal to Iq. In order to calculate self-energy diagram 
insertions to the effective Green functions in the two-loop approximation, we use the ladder 
expansion (figure [TOj) . At small values of wavevectors the bare Fq- vertex (figure [TOk) is 

= P[I{k -q) + i{k - s) - I{s) - I{q)] = 2f3Ioa\q, s), 

where 1,2; 3, 4 is the abridged notation of 2-dimensional wavevectors, which are variables of 
Fo-vertex; \k\, \q\, \s\ ^ a~^; 

Hq) = J2 ^i^^y - ^'xy) exp[-iq{Uy - I'^y)] = 2/o[cos(g^a) + cos(gya)]. 

X — f' 

The F- vertex in the ladder approximation (figure [TOb) is determined by the relationship 
F(l, 2; 3, 4) = V{k, Wi, s + q- k, u^ + uja- ^i', q, ^3, s, U4) = To{k, s + q- k; g, s) 

1 f 

-— J2 J To{k,s + q- k- q', s + q- q')G-+{q', u^^^)G-+{s + q- q', + ^4 - c^L'^) 



,(9) 
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xr('?', s + q- q', u-i + Ui- cu^); g, W3, s, ^4) £q\ 

where 

25(p) 



f3h{ijj{q) — iuj„ 



is the effective Green function determined by the P-matrix (fT2!) . (X'(g) is the frequency of 
spin excitations in monolayer (fT9|) . S'b is the volume of the 2-dimensional first Brillouin zone. 
The coefficient l/8B'^{p) is due to the fact that the substitution of the bare Green function 
to effective ones in diagrams are performed inside blocks. The self-energy diagram insertion 
(figure [TUb ) is given by 



ro(g, s + k - g; s, - + g, - ^ + 0;^^'^) 



X G^^ik, a;r))r(.-', ^, g-, o;^^), - q, u;f + cuf^) - ul^^) d'k <fs. (39) 

The damping of spin wave excitations is expressed by the imaginary part of the self-energy 
n(g-',a;(^)^ 



(40) 

iuJm -^u}+i£ sign uj 



uj{q) (3oj 

Taking into account the self-energy Ii[q,uj'^) in the Born approximation, namely, substituting 
r — )■ To in relation (!39l) . integrating over k^ s and summing over the frequency variables w^^^ 
and ujf'' in equation (HUi) . at huj{q) < kT we obtain 

C{qa)\kTf 



167r52(p)Joe(o)' 

where C = 1.12, /c is the Boltzmann constant, e^^-* = h'y(H + Z/^™)) is the Zeeman energy. 
In order to evaluate the damping of spin waves, we calculate A(g) for spin waves with the 
wavelength A = 5 fim propagating in the monolayer film with the lattice constant a = 0.4 nm 
and with the exchange interaction between neighboring spins Iq = 0.085 eV, B{p) = 1/2 at 
T = 300 K. Then, taking into account that g = 27r/A, for e^^^/h = 10 GHz we obtain A(g) = 
4.28-10"^. Thus, one can see that the damping of spin waves of small wavevectors is low. 

6 Spin-wave devices on the base of nanosized magnetic 
films 

According to the above-mentioned, spin excitations in nanosized films have low damping. This 
property can be used in spin-wave devices. We consider spin-wave filters on the base of nano- 
sized magnetic films and field-effect transistors with magnetic films under the gate contact. 
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Figure 10: (a) Bare Fq- vertex, (b) Ladder approximation, (c) Self-energy diagram insertion. 
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6.1 Spin- wave filters 



Using thin magnetic films, we can construct spin-wave devices of small sizes and can integrate 
them to semiconductor chips. Figure [TTk presents tunable spin-wave filter on a piezoelectric 
substrate. Microwave frequency current flowing in the microstrip antenna of the width w 
generates spin waves in the magnetic film. Wavevector of spin waves is determined by the 
antenna width and is in the range [0, 271/11)]. When spin waves come up to the second antenna, 
the magnetic field of spin waves induces a current of the same frequency in this antenna. The 
waveguide impedance of the filter depends on the antenna width, the width of microstrips, 
the dielectric constant of the film, and the thickness of the film between microstrips and the 
upper metal contact on the piezoelectric substrate. In order to remove a lattice mismatch and to 
reach desirable impedance, buffer layers between the magnetic film and the metal contact can be 
used. Tunability of the filter is provided by lattice variations of the piezoelectric substrate. The 
applied voltage U varies the lattice constant of the substrate and, as a result of this variation, 
varies the lattice constant of the magnetic film. Compression and expansion of the lattice 
lead to the stress anisotropy H^"''> in the magnetic film [38j- In this case, we must substitute 
H ^ H + in the spin propagators ([7]), in the Green functions (ITT!) , in the frequency fl{q) 
in ( !T9|) . ( 120|) . and in the Landau-Lifshitz equations ( |28|) . (!29|) . The frequency u of spin waves 
is varied. 

Narrow-band filters can be constructed on the base of periodic antennae (FigurefTTb). These 
antenna structure generate spin waves with the wavevector q = 27r//, where / is the period 
of the generating antenna. Filters with periodic antenna structures have more selectivity in 
comparison with filters with single antennae. The bandwidth of the filter is given by 

2n 



where u{q) is dispersion relation (fT9i) . (!20l) . or ( |33i) . go = 2it/1, Ni, N2 are numbers of periods 
of generating and receiving antennae, respectively. Thin magnetic film used in the above- 
mentioned filters must be dielectric and can be garnet, spinel or hexaferrite films. These films 
can be produced by the laser deposition method or by ion-beam sputtering with following 
annealing. In [SH] band pass spin-wave filters at 5 — 7 GHz have been fabricated on the base 
of submicron thick YIG films produced by the laser deposition on Gd3Ga50i2 substrates. In 
the case, when spin-wave filters are integrated on semiconductor chips, the semiconductor 
must endure the annealing procedure without any changes for the worse in the semiconductor 
structure. For this purpose. Si and GaN can be used. 

6.2 Field-effect transistors with nanosized magnetic films 

As we can see from relations (120|) . (pTj) and from figures HJ [5], spin- wave resonance peaks in thin 
magnetic films have high frequencies and low numbers. For the exchange interaction between 
layers Id = 0.05-0.1 eV (garnet films, lithium ferrospinel and hexaferrite films [3SII1D]) spin-wave 
resonance frequencies are in Gigahertz and Terahertz frequency bands. Distance between spin- 
wave resonance peaks is greater in thin films than in thick ones. Since in thin films spin-wave 
resonances fl2Tl) at the given high frequency have lower numbers n and have lower numbers of 
half-periods of waves on film thickness in comparison with resonances at the same frequency in 
thick magnetic films, high-frequency resonances in thin films can be exited more easily. This 
can be used in field-effect transistor (FET) structures. In figure [12] thin dielectric magnetic 
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Figure 11: (a). Spin- wave filter on the base of nanosized magnetic film on a piezoelectric 
substrate, (b). Spin- wave filter with nanosized magnetic film on a semiconductor substrate. 
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Figure 12: Field-effect transistor with nanosized magnetic film under the gate contact. 



film is placed under the gate electrode. Applied voltage of the microwave frequency on the 
gate generates electromagnetic field, which induces spin-wave resonances. In order to increase 
the magnetic component of the electromagnetic field and to enhance amplitude of resonances, 
multiferroic layer can be placed between the gate and the magnetic film. The magnetic field of 
the induced spin-wave resonances interacts with spins of electrons propagating in the channel 
and modulates the current of the transistor. Consequently, existence of spin-wave resonances 
in the gate circuit can be regarded as a filter. Choosing the certain spin-wave resonance peak, 
we can construct FET structure operating in a desirable range of Gigahertz and Terahertz 
frequency bands. 



7 Conclusions 

The results of the investigations can be summarized as follows. 

(1) We have studied spin excitations in nanosized magnetic films in the Heisenberg model with 
magnetic dipole and exchange interactions by the spin operator diagram technique. Dispersion 
relations of spin waves in two-dimensional magnetic monolayers and in two-layer magnetic films 
and the spin-wave resonance spectrum in A^-layer structures are found. 

(2) Generalized Landau-Lifshitz equations for thick magnetic films which are derived from first 
principles, have the integral (pseudodifferential) form, but not differential one with respect to 
spatial variables. Spin excitations are determined by simultaneous solution of the Landau- 
Lifshitz equations and the equation for the magnetostatic potential. The use of exchange 
boundary conditions for solvability of the Landau-Lifshitz equations is incorrect. 

(3) The magnetic dipole interaction makes a major contribution to the relaxation of long- 
wavelength spin waves in thick magnetic films. The spin-wave damping is determined by 
diagrams in the one-loop approximation, which correspond to three-spin-wave processes. The 
three-spin-wave processes are accompanied by transitions between thermal excited spin-wave 
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modes. The damping increases directly proportionally to the temperature. 

(4) Thin films have a region of low relaxation of long-wavelength spin waves. In thin magnetic 
films the energy of these waves is less than energy gaps between spin-wave modes, therefore, 
three-spin-wave processes are forbidden, four-spin-wave processes take place and, as a result of 
this, the exchange interaction makes a major contribution to the relaxation. It is found that 
the damping of spin waves propagating in a magnetic monolayer has the form of the quadratic 
dependence on the temperature and is very low for spin waves with small wavevectors. 

(5) Nanosized magnetic films can be used in spin- wave devices. Low damping of long- wavelength 
spin waves gives us opportunity to construct tunable narrow-band spin-wave filters with high 
quality at the microwave frequency range. Spin-wave resonances of nanosized magnetic films 
under gate electrodes in field-effect transistor (FET) structures can be used to construct FET 
structures operating in Gigahertz and Terahertz frequency bands. 
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